Ultra-conserved sequences in the genomes of highly diverse Anopheles mosquitoes, with implications for malaria vector control.

DNA sequences that are exactly conserved over long evolutionary time scales have been observed in a variety of taxa. Such sequences are likely under strong functional constraint and they have been useful in the field of comparative genomics for identifying genome regions with regulatory function. A potential new application for these ultra-conserved elements (UCEs) has emerged in the development of gene drives to control mosquito populations. Many gene drives work by recognizing and inserting at a specific target sequence in the genome, often imposing a reproductive load as a consequence. They can therefore select for target sequence variants that provide resistance to the drive. Focusing on highly conserved, highly constrained sequences lowers the probability that variant, gene drive-resistant alleles can be tolerated. Here, we search for conserved sequences of 18 bp and over in an alignment of 21 Anopheles genomes, spanning an evolutionary timescale of 100 million years, and characterize the resulting sequences according to their location and function. Over 8000 UCEs were found across the alignment, with a maximum length of 164 bp. Length-corrected gene ontology analysis revealed that genes containing Anopheles UCEs were over-represented in categories with structural or nucleotide-binding functions. Known insect transcription factor binding sites were found in 48% of intergenic Anopheles UCEs. When we looked at the genome sequences of 1142 wild-caught mosquitoes, we found that 15% of the Anopheles UCEs contained no polymorphisms. Our list of Anopheles UCEs should provide a valuable starting point for the selection and testing of new targets for gene-drive modification in the mosquitoes that transmit malaria.

A CRISPR-Cas9 sex-ratio distortion system for genetic control.

Genetic control aims to reduce the ability of insect pest populations to cause harm via the release of modified insects. One strategy is to bias the reproductive sex ratio towards males so that a population decreases in size or is eliminated altogether due to a lack of females. We have shown previously that sex ratio distortion can be generated synthetically in the main human malaria vector Anopheles gambiae, by selectively destroying the X-chromosome during spermatogenesis, through the activity of a naturally-occurring endonuclease that targets a repetitive rDNA sequence highly-conserved in a wide range of organisms. Here we describe a CRISPR-Cas9 sex distortion system that targets ribosomal sequences restricted to the member species of the Anopheles gambiae complex. Expression of Cas9 during spermatogenesis resulted in RNA-guided shredding of the X-chromosome during male meiosis and produced extreme male bias among progeny in the absence of any significant reduction in fertility. The flexibility of CRISPR-Cas9 combined with the availability of genomic data for a range of insects renders this strategy broadly applicable for the species-specific control of any pest or vector species with an XY sex-determination system by targeting sequences exclusive to the female sex chromosome.

Genomic signatures of population decline in the malaria mosquito Anopheles gambiae.

BACKGROUND: Population genomic features such as nucleotide diversity and linkage disequilibrium are expected to be strongly shaped by changes in population size, and might therefore be useful for monitoring the success of a control campaign. In the Kilifi district of Kenya, there has been a marked decline in the abundance of the malaria vector Anopheles gambiae subsequent to the rollout of insecticide-treated bed nets. METHODS: To investigate whether this decline left a detectable population genomic signature, simulations were performed to compare the effect of population crashes on nucleotide diversity, Tajima's D, and linkage disequilibrium (as measured by the population recombination parameter ρ). Linkage disequilibrium and ρ were estimated for An. gambiae from Kilifi, and compared them to values for Anopheles arabiensis and Anopheles merus at the same location, and for An. gambiae in a location 200 km from Kilifi. RESULTS: In the first simulations ρ changed more rapidly after a population crash than the other statistics, and therefore is a more sensitive indicator of recent population decline. In the empirical data, linkage disequilibrium extends 100-1000 times further, and ρ is 100-1000 times smaller, for the Kilifi population of An. gambiae than for any of the other populations. There were also significant runs of homozygosity in many of the individual An. gambiae mosquitoes from Kilifi. CONCLUSIONS: These results support the hypothesis that the recent decline in An. gambiae was driven by the rollout of bed nets. Measuring population genomic parameters in a small sample of individuals before, during and after vector or pest control may be a valuable method of tracking the effectiveness of interventions.

Partitioning the contributions of alternative malaria vector species.

BACKGROUND: In many locations malaria is transmitted by more than one vector species. Some vector control interventions, in particular those using genetic approaches, are likely to be targeted against a single species or species complex, at least initially, and it would therefore be useful to be able to predict the epidemiological impact of controlling a single species when multiple vector species are present. METHODS: To address this issue, the classical Ross-McDonald model of malaria epidemiology is expanded to account for multiple vector species, giving expressions for the equilibrium prevalence, sporozoite rates and reproductive number. These allow one to predict when control of just one vector species will lead to elimination of the disease. Application of the model is illustrated using published data from a particularly extensive entomological and epidemiological survey before the rollout of bed nets in eastern Kenya, where Anopheles gambiae s.l. and An. funestus were vectors. RESULTS: Meta-analysis indicates that sporozoite rates were 38 % higher in An. gambiae s.l. than in An. funestus, and, according to the model, this difference could be due to An. gambiae s.l. having a higher frequency of feeding on humans, a higher human-to-mosquito transmission rate, a lower adult mortality rate, and/or a shorter incubation period. Further calculations suggest that An. gambiae s.l. would have been sufficient to maintain transmission by itself throughout the region, whereas An. funestus would not have been able to support transmission by itself in Malindi District. CONCLUSIONS: Partitioning the contributions of different vector species may allow us to predict whether malaria will persist after targeted vector control.

Genomic analyses of three malaria vectors reveals extensive shared polymorphism but contrasting population histories.

Anopheles gambiae s.l. are important malaria vectors, but little is known about their genomic variation in the wild. Here, we present inter- and intraspecies analysis of genome-wide RADseq data, in three Anopheles gambiae s.l. species collected from East Africa. The mosquitoes fall into three genotypic clusters representing described species (A. gambiae, A. arabiensis, and A. merus) with no evidence of cryptic breeding units. Anopheles merus is the most divergent of the three species, supporting a recent new phylogeny based on chromosomal inversions. Even though the species clusters are well separated, there is extensive shared polymorphism, particularly between A. gambiae and A. arabiensis. Divergence between A. gambiae and A. arabiensis does not vary across the autosomes but is higher in X-linked inversions than elsewhere on X or on the autosomes, consistent with the suggestion that this inversion (or a gene within it) is important in reproductive isolation between the species. The 2La/2L+(a) inversion shows no more evidence of introgression between A. gambiae and A. arabiensis than the rest of the autosomes. Population differentiation within A. gambiae and A. arabiensis is weak over approximately 190-270 km, implying no strong barriers to dispersal. Analysis of Tajima's D and the allele frequency spectrum is consistent with modest population increases in A. arabiensis and A. merus, but a more complex demographic history of expansion followed by contraction in A. gambiae. Although they are less than 200 km apart, the two A. gambiae populations show evidence of different demographic histories.

Genetic population structure of the malaria vector Anopheles baimaii in north-east India using mitochondrial DNA.

BACKGROUND: Anopheles baimaii is a primary vector of human malaria in the forest settings of Southeast Asia including the north-eastern region of India. Here, the genetic population structure and the basic population genetic parameters of An. baimaii in north-east India were estimated using DNA sequences of the mitochondrial cytochrome oxidase sub unit II (COII) gene. METHODS: Anopheles baimaii were collected from 26 geo-referenced locations across the seven north-east Indian states and the COII gene was sequenced from 176 individuals across these sites. Fifty-seven COII sequences of An. baimaii from six locations in Bangladesh, Myanmar and Thailand from a previous study were added to this dataset. Altogether, 233 sequences were grouped into eight population groups, to facilitate analyses of genetic diversity, population structure and population history. RESULTS: A star-shaped median joining haplotype network, unimodal mismatch distribution and significantly negative neutrality tests indicated population expansion in An. baimaii with the start of expansion estimated to be ~0.243 million years before present (MYBP) in north-east India. The populations of An. baimaii from north-east India had the highest haplotype and nucleotide diversity with all other populations having a subset of this diversity, likely as the result of range expansion from north-east India. The north-east Indian populations were genetically distinct from those in Bangladesh, Myanmar and Thailand, indicating that mountains, such as the Arakan mountain range between north-east India and Myanmar, are a significant barrier to gene flow. Within north-east India, there was no genetic differentiation among populations with the exception of the Central 2 population in the Barail hills area that was significantly differentiated from other populations. CONCLUSIONS: The high genetic distinctiveness of the Central 2 population in the Barail hills area of the north-east India should be confirmed and its epidemiological significance further investigated. The lack of genetic population structure in the other north-east Indian populations likely reflects large population sizes of An. baimaii that, historically, were able to disperse through continuous forest habitats in the north-east India. Additional markers and analytical approaches are required to determine if recent deforestation is now preventing ongoing gene flow. Until such information is acquired, An. baimaii in north-east India should be treated as a single unit for the implementation of vector control measures.

Comparative phylogeography reveals a shared impact of pleistocene environmental change in shaping genetic diversity within nine Anopheles mosquito species across the Indo-Burma biodiversity hotspot.

South-East Asia is one of the world's richest regions in terms of biodiversity. An understanding of the distribution of diversity and the factors shaping it is lacking, yet essential for identifying conservation priorities for the region's highly threatened biodiversity. Here, we take a large-scale comparative approach, combining data from nine forest-associated Anopheles mosquito species and using statistical phylogeographical methods to disentangle the effects of environmental history, species-specific ecology and random coalescent effects. Spatially explicit modelling of Pleistocene demographic history supports a common influence of environmental events in shaping the genetic diversity of all species examined, despite differences in species' mtDNA gene trees. Populations were periodically restricted to allopatric northeastern and northwestern refugia, most likely due to Pleistocene forest fragmentation. Subsequent southwards post-glacial recolonization is supported by a north-south gradient of decreasing genetic diversity. Repeated allopatric fragmentation and recolonization have led to the formation of deeply divergent geographical lineages within four species and a suture zone where these intraspecific lineages meet along the Thai-Myanmar border. A common environmental influence for this divergence was further indicated by strong support for simultaneous divergence within the same four species, dating to approximately 900 thousand years ago (kya). Differences in the geographical structuring of genetic diversity between species are probably the result of varying species' biology. The findings have important implications for conservation planning; if the refugial regions and suture zone identified here are shared by other forest taxa, the unique and high levels of genetic diversity they house will make these areas conservation priorities.

Molecular phylogenetics and biogeography of the Neocellia Series of Anopheles mosquitoes in the Oriental Region.

Molecular studies of population divergence and speciation across the Oriental Region are sparse, despite the region's high biodiversity and extensive Pliocene and Pleistocene environmental change. A molecular phylogenetic study of the Neocellia Series of Anopheles mosquitoes was undertaken to identify patterns of diversification across the Oriental Region and to infer the role of Pleistocene and Pliocene climatic change. A robust phylogeny was constructed using CO2 and ND5 mitochondrial genes and ITS2 and D3 nuclear ribosomal markers. Bayesian analysis of mitochondrial genes was used to date divergence events. The repeated contraction and expansion of forest habitat resulting from Pleistocene climatic fluctuations appears to have had a substantial impact on intraspecific diversification, but has not driven speciation within this group. Primarily early to mid Pliocene speciation was detected within the Annularis Group, whereas speciation within the Maculatus and Jamesii Groups occurred during the mid and late Pliocene. Both allopatric divergence driven by late Pliocene environmental changes and ecological adaptation, involving altitudinal replacement and seasonality, are likely to have influenced speciation in the Maculatus Group.

Crossing experiments supporting the specific status of Anopheles maculatus chromosomal form K.

There are 3 recognized chromosomal forms (B, E, K) in the taxon of Anopheles maculatus, 1 of the 8 members of the Anopheles maculatus group. Previous studies suggested that forms B and E are cytotypes of the species, but genetic characteristics of form K are unknown. The present study used crossing experiments, and polytene chromosomes of the ovarian nurse cell in F1 hybrids to show that form K is genetically distinct from forms B and E. In addition, postzygotic genetic incompatibility between form K and An. sawadwongporni, An. dravidicus, and An. pseudowillmori are demonstrated. In all crosses, hybrid males were sterile, with atrophied testes and accessory glands, or partially sterile with abnormal spermatozoa. The hybrid females showed varying degrees of atrophied ovaries. The ovarian nurse cell polytene chromosomes of the F1 hybrid females displayed approximately 70% to almost complete asynapsis. The results provide clear evidence that form K should no longer be regarded as a cytotype of An. maculatus. The present study supports previous suggestions that form K represents another species of the Maculatus group.